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FOREWORD This Indian Standard (Part 3/Set 2) was adopted by the Bureau of Indian Standards, after the draft finalized by the Electrical Insulation Systems Sectional Committee had been approved by the Electrotechnical Division Council. This standard is one in the series of standard dealing with guidelines for the evaluation of insulation systems of electrical equipment. These guidelines include method of identification and general principles of functional evaluation of insulation systems together with test procedures to evaluate the behaviour of systems under different factors of influence. The purpose of this standard is to guide the development of system test procedures and suggest pcints to be considered by the equipment technical committees in the preparation of specific instructions for the evaluation of electrical endurance of insulation systems for electrical equipment. This standard suggests techniques that may be used for data whose distribution does not correspond well with a normal (Gaussian) distribution, but may correspond to an extreme-value distribution. The other parts of the standard are as follows: IS 11182 Part 1 Part 2 Guide for the evaluation of electrical insulation system of electrical equipment: Identification, evaluation and ageing mechanism Thermal endurance test procedures and

Part 3/ Set 1 Electrical endurance test procedures, Section 1 General considerations evaluation procedures based on normal distribution Part 5 Part 6 Mechanical endurance functional tests Performance evaluation based on service experience and functional tests

Part 7/ Set 1 Multi factor functional testing, Section 1 Test procedures In the preparation of this standard, assistance has been derived fron IEC Publication 727-2 (1993) `Evaluation of electrical endurance of electrical insulation system: Part 2 Evaluation procedures based on extreme-value distribution', issued by the International Electrotechnical Commission (IEC). For the purpose of deciding whether a particular requirement of this standard is complied with, the final value, observed or calculated, expressing the result of a test or analysis, shall be rounded off in accordance with IS 2 : 1960 `Rules for rounding off numerical values ( revised )' . The number of significant places retained in the rounded off value should be the same as that of the specified value in this standard.
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Indian Standard

GUIDEFORTHEEVALUATIONOFINSULATION SYSTEMSOFELECTRICALEQUIPMENT
PART 3 ELECTRICAL ENDURANCETEST PROCEDURES Section 2 Evaluation Procedures Based on ExtremsValue Distributions IS No.

1 SCOPE 1.1 This standard describes statistical procedures to analyze data on time-to-breakdown at constant voltage stress (or breakdown voltages under long-term increasing-voltage stress) of groups of individual specimens of electrical insulation systems, or models representinginst&ingsystems. Numericalexamplesare also included.

me
evaluation Performance based on service experience and functional tests MuIti factor functional Section 1 Test testing, procedures

(Part6) : 1986

(Part7/Sec

1) : 1986

1.2 This standard assumes that the electrical factors of influence is the dominating ageing factor. If it is known or understood that more than one factor of influence is significant, reference should be made to IS 11182 (Part 7/ Set 1) : 1986 `Guide for evaluation of insulation systems of electrical equipment: Part 7 Multi factor functional testing, Section 1 Test procedures'. 1.3 This standard does not apply to data derived from short-time electrical breakdown tests except those used as diagnostic factors from other long-term factors of intluence. When used in this way for diagnosis of factors of influence other than electrical, reference should be made to other appropriate Indiau Standards given in 2 of this standard. 2 REFERENCES The following Indian Standards are the necessary adjuncts to this standard :
IS No. Title

3 PROBABIUI'Y DLSTRIBUTIONS The statistical analysis of a set of data presupposes that a probability distribution function has been chosen that is assumed to represent the variation of the data under the conditions of the test. In this guide two types of extreme-value distribution are described. 3.1 Extrem~Value
Distributions

Extreme-value distributions have been found to represent a wide variety of data on systems that fail by a `weakest-link' mechanism, and they are often employed to analyze time-to-breakdown and long-time breakdown-voltage data in insulation studies. The extreme-value distributions are fundamentally ditferent from the normal ( Gaussian ) distribution. This latter distributionbest represents `average' phenomena. Techniques commonly used for estimating normal distribution quantities are not generally applicable to the extreme-value distributions. In particular, the methods to calculate the mean, standard deviation, confidence intervals and hypothesis quantities for the normal distribution should not be applied to extremevalue distributions, Two types of extreme-value distributions Weibull and the Gumbel distributions. 3.1.1 Weibull Distribution
The Weibull distribution is appropriate when the rate

2071 (Part 1) : 1993

High voltage test techniques: Part 1 General definitions and test requirements (second revision ) Guide for evaluation insulation systems electrical equipment: of of

are the

11182

(Part 1) : 1984 (Part2) : 1984

Identification, evaluation and ageing mechanism Thermal endurance procedures test

of failure varies with time. It is most often used to represent either time-to-failure results from a longtime constant voltage stress test, or breakdown voltage results form a long-time progressive voltage stress test, on solid insulations. The three-parameter Weibull distribution has the
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3/Set 2 ) : 1996 divergence of the true graph from a straight line will increase as the difference between the true value and the assumed value of the location parameter (r) increases. The two-parameter following equation:
F(x) =.I--exp

following equation:

F(x)
where
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Weibull distribution

has the

scale parameter and is positive; shape parameter and is positive; location parameter; random variable, usually the time to breakdown or the breakdown voltage; and of 5x. failure at time
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3.1.2 Gumbel Distribution
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F(x) = probability (or voltage)

The probability of failure F(x) is zero at x < y The probability of failure rises continuously as x increases. As the time or voltage increases to infinity, the probability of failure approaches certainly, that is F(x=oo)=l. The scale parameter (a) represents the time (or voltage) required for (l- e-l) or 63.2 percent of the tested specimens to fail. It is analogous to u , the mean of the normal distribution*. The unit of a is the same as that of x, that is, hours or volts. The shape parameter (p) is a measure of dispersion of the log failure times or log voltages. The larger p is, the smaller is the range of log times or log breakdown voltages. The shape parameter, which is dimensionless, is analogous to l/a, where cr is the standard deviation of the normal distribution*. When the shape parameter (p) is less than 1, the failure rate decreases with time; when greater than 1, the faihrre rate increases with time. When the shape parameter (p) is equal to 1, the Weibull distribution is equivalent to the exponential distribution and the failure rate is independent of time. The location parameter (y) represents a value below which the probability of occurrence of x is zero. It represents a threshold value, below which a certain physical phenomena (failure) does not occur. The unit of y is the same as that of x. The estimation of three Weibull parameters requires the solution of non-linear equations, and is not possible by a two-dimensional graphical solution. The two-parameter Weibull distribution is a particular case of the three-parameter distribution, when the locationparameter (r) is zero. This latter implies that the phenomenon being studied has a threshold value equal to zero. One a two-dimensional graphical representation, the
*Note that a mean and a standard deviation can be defined for the Weibull distribution. These quantities however, are seldom used.

The Gumbel distribution is most often used to represent the breakdown voltages of liquid insulation and of compressed gas insulations with slightly nonuniform fields. The cumulative Gumbel distribution function for the population fraction below y is:

where u = location parameter and may have any
VdUe;

b = Y= GO

scale parameter

and is positive;

random variable, usually the breakdown voltage; and = probability of failure at voltage (or time) less than or equal toy.

The Gumbel distribution is unsymmetrical and has a physically impossible finite probability of breakdown for y < 0, but if u >> b, this probability will be negligibly small. This distribution is also called the smallest extreme-value (that is, weakest-link) distribution. The unitof uandbisthesameasthat of y. The Gumbel distribution is closely related to the Weibull distribution. That is, ifx has a Weibull distribution, then y = In (x) has a Gumbel distribution where: u = Ina ......(4)

Estimation techniques employed for one distribution ( Gumbel or Weibull ) apply to the other if the transformation equations (4) and (5) are utilized. 4 TREATMENTOFTFSTDATA 4.1 Censored Data Censored data occur when n specimens are started on test together and the times to breakdown of only r (< n) are observed. Censoring is encountered mainly with constant voltage stress tests, where the data are 2

IS 11182 ( Part 3/Set 2 ) : 1996

analyzed or the test is terminated before all the specimens fail. Censoring can also occur with progressive-stress tests where f&hovers or spurious breakdowns occur on specimens which survive to high stresses. Since censoring can occur by plan or by accident in many insulation tests, it should be taken into account in the data analysis. There are at least two types of censoring: Typel: Occurs when the test is terminated (or the data are analysed) after a certain period of time (or given maximum test voltage ). Suppose x, is the time the test is stopped and xr is the time of therth and last observed failure (x, > x, ), then (n -r) unfailed specimens survived beyond time x, Occurs when the test is terminated immediately afler the rth failure, that is, x,=x,.

voltages from smallest to largest. The cumulative probability of the ith smallest value (xi) is approximated by the formula :
F(Xi)=n+l .....

where n is the total number of test specimens. The Weibull example data in Table 1 are plotted in Fig. 1A. Data from unfailed specimens should not be plotted, although they should be included in n. Specimens that have failed by mechanisms that are clearly spurious and invalid, should be disqualified, and not counted inn. If the plotted data follow a straight line, then it may be reasonable to assume that the results are adequately represented by the Weibull distribution. Some random deviations from a straight line would normally be expected. If, however, there is a consistent departure from a straight line (say, curvature), then another distribution may fit the data better. Figure 2 shows Weibull probability plots of data from other distributions. 4.2.2 Gumbel Probability Graph Paper A sample of Gumbel probability paper is shown in Fig. 3A. As with the Weibull paper, the failure times or voltages are ordered from smallest to largest, and the cumulative probability for the ith smallest value @i) is approximated by the formula:

TypeII:

There are also other kinds of censoring, such as progressive and multiple censoring. These types result, for example, when failed specimens are disqualified because of specimen failure by spurious mechanisnis. Analytical treament of such data is more diflicu1t. 4.2 Selection of a Probability Distribution Solid insulation test data are often represented by the two-parameter Weibull distribution whereas the Gumbel ( smallest extreme-value ) distribution is often employed for liquids. However, the validity of the assumed distribution should be tested since it is possible that another distribution will yield a better fit, unless prior experience or theory indicates that a particular distribution is valid. Testing the validity of the assumed distribution to best represent the test data can most easily be accomplished by the use of probability graph paper. This paper has one axis in the form of a non-linear cumulative probability scale, which is the same for both the Weibull and the Gumbel distributions. The other axis, on which are plottedthe failure times or breakdown voltage, is lograithmic for the Weibull distribution and linear for the Gumbel distribution. The axis are scaled so that data from the respective distributions will follow a straight line on the corresponding graph paper. Other methods to test the validitv of assumed distributions include statistical fitness tests, or comparison by the maximum likelihood method.
4.2.1 Weibull Probability Graph Paper

G(Yi) =

I.

EAT
~ I

.....

where n is the totoal number of test specimens. If the data plot as a straight line on the Gumbel paper, then the Gumbel distribution adequately represents the failure data. Examples of data from other distributions are plotted on Gumbel paper in Fig. 4. 4.2.3 Graph Paper Availability The graph papers described in 4.2.1 and 4.2.2 are generally commercially available. However, if double natural-logarithmic transformations are made, the Weibull function becomes linear with respect to In x, and the data may be plotted on standard arithmetic paper. Appropriate data will form straight lines. 4.3 Estimation of Distribution Parameters The distribution parameters may be estimated graphically or by calculation. The calculated estimates are more objective than the graphical estimates and are preferred. 4.3.1 Graphical Estimation of Parameters The probability graph papers described in 4.2 can be used to obtain approximate estimates of the parameters, from the slope and intercept of the straight line fitted to the data points `by eye'. 3

A sample of Weibull probability paper is shown in Fig. 1A. To use this paper, order the failure times or

/

Origin 99-9
98 90

70
50 30

20 IO 5*0

2-o

1'0 0'5

0'2 0'1 30 Time (hours)

100

1 000

FIG. 1A WEIBULL PLOTOFFAILURE DATAWITH LINE FIT `BY EYE'

999
98 90 70 50 30

_................................................................................................... i.................................................

.....................................................

.

.....................................................................

....................................................................................................

2'0 1'0

.....................

.

.........

..,

.....................................................................

0'5

_.

..................................................................................................

i .................................................

.....................................................

i

.....................................................................

0'2 0'1 1

I
10 Time (hours) FIG. 30

I
100 500

1B WEIBULLPLOT OF90 PERCENT CONFIDENCE BOUNDSFORPERCENTILES AND MAXIMUM LIKELIHOOD FITTED LINE

999 98 90 70 50 30 20

l-

/

Sparameter

Weihull

/

lormal

/

/ 10 Gumbel 5'0 /

-

2'Cb -

/

l'( I/ P` .S-

/
O'I2-

/

0' 1-

Time (hours)

FIG. 2 PLOT OF OTHERDISTRIBUTIONS ON~-PARAMETER WEIBULL PAPER

IS 11182 (Part 3/Set 2) : 1996 Confidence intervals for parameters cannot be estimated graphically.
4.3.1.1 Graphical estimatesfor the Weibull distribution

However, for large sample sizes (for example n > 50) of uncensored data, transformed as described in 4.2.3, the errors inherent in linear regression methods, such as the least-squares method, are not usually apparent. 4.3.2.1 Calculated estimatesfor the Weibull distribution For the Weibull distribution, the maximum likelihood estimates of & and b require the iterativesolutions of: f(p)= and 2fC=Oand .....(n)

Plot the test data on Weibull probability paper, as described in 4.2.1. Fit a striaght line `by eye' to the data points. The estimate for the scale parameter ((r), denoted by (8) is the time (or voltage) corresponding to F(x) = 63.2 perce$. An estimate for the shape parameters (p) denoted by (p), is calculated from: b= In (ln L1) - In (In b) ......(S) In ( z where L, = 1 1- F(x,) ......(lO) 1F(xz)

> where ......(9) A, 7 $x? [ln .Qk-`) + (n-r) x,*[ln x.]`"`)

......(12)

...... (13)

+'

for k =I, 2,3 (A3 is required for the iterative solution).

In ( ) represents the natural logarithm and [xl, F(x,)J, [x2, F(xx)] correspond to two points on the fitted line. Most commercial Weibull probability papers contain a special scale which permits the rapid estimationof p. For the data from Table 1 plotted in Fig. 1A, & =130 h. Usingthetwopoint# (160,0.7)and (30,0.16) obtained from the fitted line, p= 1.15 using equation (8). The line labelled A-A, drawn through a point labelled `origin' in Fig. 1A, is parallel to the fitted line and gives an estimate of 1.15 for p. 4.3.1.2 Graphical estimatesfor the Gumbel distribution Plot the test data on Gumbel probability paper, as described in 4.2.2. Fit a straight line `by eye' to the data points. The estimate of the location parameter (u), denoted by (u^), is the voltage corresponding to G(Y)= 63.2 percent. An estimate of the scale parameter (b), denoted by (8), is the difference between the 63 percent and 3 1 percent points. For the data from Table 2 plotted in Fig. 3A, u^ = 5.84 kV Using the two percentage points, from the fitted line, 8 = 5.84 - 5.43 = 0.41 kV. 4.3.2 Calculated Estimation of Parameters The most widely used parameter estimates are the maximum likelihood estimates, which are valid for both Type I and Type II censored data. They are asymptotically unbiased and their variances are not larger than the variances of any other estimates, for example those derived from the least-squares method.

where x, = running time or highest test voltage of the experimerd, = total number of specimens tested, and n r = number of specimens that failed. For Qpe II[Icensoring, xS=x,. A short BASIC programme solving these equations on a personal computer is shown in Annex A. Most computer programs require an initial guess for p. This is most easily obtained from a Weibull probability plot. In lieu of probability plot, an initial p of 10 may be applicable for progressive-stress tests. For constantstress tests, an initial p of 2 may be applicable. The maximum likelihood estimates for the data in Table 1 are cr^ =115 hand b=1.5. The line on Weibull probability paper which fits these data is, shown in Fig. 1B. The line is plotted by substituting two values of x, along with 8 and b into equation (2) and calculating the corresponding values of F(x). 4.3.2.2 Calculated estimatesfor the Gumbel distribution A separate set of equations to yield the maximum likelihood estimates could be generated for the Gumble distribution . However, it is simpler to employ the
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,IS 11182 ( Part 3/See 2 ) : 1996 transformations in equations (4) and (5) with a computer program that fits the Weibull distribution. In particular, the Gumbel test data wi) are first transformed into Weibull data with:
Xi= eXp(_Yi)

Many tables similar to Student's t table have been calculated for the extreme-value distributions. For those without access to large computers, approximate confidence intervals based on these tables can be determined. 4.4.1 Approximate Parameters
Confidence Intervals for WeibuN

......(15)

A first estimate for p = i can be obtained from a Gumbel probability plot. Alternatively, assume an initial b of 0.1 (that is, p =lO). The programme output of (r,and 8 are then: u^=ln^a ...... (16)

There are various methods and tables for calculating approximate contidence intervals for 01 and l3.The most complete tables are based on best linear invariant estimates for the parameters, rather than the maximum likelihood estimates. The confidence intervals tables for a and 8, which cover all possible censoring fractions are very extensive. For simplicity these tables are represented as curves in Fig. 5 and 6. These curves are approximate, valid only for Type II censoring, and are for experiments with up to 25 specimens tested. Type I censored data can be analyzed with some bias in the estimates by setting x, = x,. Figure 5 and 6 give 90 percent confidence intervals only For more accurate intervals, other sample sizes, etc., computer-based methods will be necessary. Figure 5 is used to compute the 90 percent confidence intervals for the shape paramenter p: Pl=K^P a=w"^P ...... (18) ...... (19)

...... (17)

For the Gumbel data in Table 2, u^= 5.73 kV and b = 0.26 kV With two values of y and the parameter estimates, corresponding values of G(y) can be calculated. These two points will define the straight line on Gumbel paper (Fig. 3B.) 4.4 Confidence Intervals for Parameters
If the same experiment involving the testing of many

specimens is performed a number of times, the values of the parameter estimates ( a^ b br a^ & ) from each experiment differ. This variation in estimates results from the statistical nature of insulation breakdown. Therefore, any parameter estimate differs from the `true' parameter value of the total population which could be obtained from an experiment involving an itiitely large number of specimens. Hence, it is common to give with each parameter estimate a `confidence interval', which encloses the true parameter value with high probability. The confidence interval objectively quantifies the uncertainty in the estimate of the parameter. In general, the more specimens tested, the narrower will be the confidence interval. Ifan experiment is poorly performed, for instance if the applied voltage is not held constant in a constant stress test, the contidence intervals are inaccura te. Confidence intervals are valid only for identically tested specimens. There are various confidence intervals for Weibull and Gumbel parameters and percentiles. The recommended intervals employ maximum likelihood parameter estimates. This method, which is vahd for all sample sizes and Type II censoring, requires an extensive computer programmes which is not generally available. Many available computer programmes give parameter estimates and approximate contidence intervals. However, their validity should be checked by comparison to reputable main-frame programmes.

where pi and & are the lower and upper bounds, respectively, for the interval. As seen in Fig. 5, WI and W, are mainly functions of r, the number of failures. For the Weibull data in Table 1, WI= 0.42 and W,= 1.46 and, since 8 1.5 (see 4.3.2.1), the 90 percent confidence limitsfor ~are~=0.42Xl.5=0.63and~=1.46~1.5 =2.2. Fig. 6 is used to calculate the 90 percent confidence intervals for the scale parameter (0~):
a,

z = ii exp 4

[ B

a" = &exp + [ B

Z

I I

......(20)

.......(21)

where aI and 041 refer to the lower and upper confidence bounds, respectively. The factor 2, is primarily a function of n, the number of specimens put on test. Z,, however, is a function of both r, the number of specimens which actually failed, and n. For the Weibull data in Table I, Z,= 0.75 and Z,=l .O.Therefore, the 90 percent amfldence limitsfor a are ai =69 hand o,,=220h.
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IS 11182 (Part J/Set 2) : 1996 4.4.2 Approximate Parameters Confidence Intervalsfor Gumbel

is:
1

The same curves used for the Weibull parameter intervals can be used for the Gumbel parameters. As before, the curves in Fig. 5 and Fig. 6 yield approximate 90 percent confidence intervals for Type II censored data. The approximate 90 percent confidence intervals for b has lower and upper limits: I\
b,

:r=;[-ln(

l-p)]

$

......(26)

where & and p^ are thg maximum-hklihood parameter estimates. The confidence interval for xPis: xl@)=&exp

I7
xU@)=Gexp

=b
W"

......(22)

bl

=b
Wl

......(23)

VI 1 VU 3 11

......(27)

.. . ..(28)

Note that the W, factor from Fig. 6 is for the lower bound and WI is for the upper bound. For the Gumbel datainTable2, WU=1.44andW1=0.45;thusb,=0.18kV and b,= 0.58 kV are the confidence limits for b. The approximate 90 percent confidence intervals for u has lower and upper limits: 111 =c-$2, Zf,=^u+22, ......(24) ......(25)

where xl (._, a) and x,, (D) refer to the lower and upper bound for the (100 b)t?h'percentile, respectively. The V factors are obtained from Fig. 7A and 7B for the first, fifth and tenth percentiles @ = 0.01, 0.05 and 0.10 respectively). V, is primarily only a function of r, whereas I',,is primarily a function of n. For the Weibull data in Table 1, for the first percentile, I',= 10.4 and VU= 3.1. Thus the 90 percent confidence limitsarex1(0.01)=0.11handxU(0.01)=14h.Thebest estimate of the first percentile is x^ r = 5.3 h. The confidence limits for the percentiles, together with the confidence interval for a can be usefully displayed on Weibull probability paper. For the upper bound, plot the calculated upper limits (x values) corresponding to the first, fifth, tenth and 63.2 ( a ) percentiles on the graph paper. Join these four points with a smooth line. Similarly, draw a line for the lower contidence limits. As shown in Fig. lB, these confidence limits will bound the `best' line calculated from the maximum likelihood parameters. On repeated tests, the best lines should predominantly fall within the confidence bounds. The greater the number of samples tested, the narrower will be the gap between the upper and lower bounds. 451.2
Confidence intervalsfor Gumbelpercentiles

The 2 factors are obtained from Fig. 6A and 6B. For the GumbeldatainTable2,Z,=0.7andZU=1.12.Thusthe 90 percent confidence intervals for u has limits ui= 5.57 kVanduU=6.04kV. 4.5 Percentiles In addition to estimating the parameters of a distribution, in many cases it is of interest to calculate the time-to-breakdown or the breakdown voltage at low probabilities of failure. The values of x and y at a particular probability of failure are referred to as percentiles ifthe probability is expressed in per cent. 4.51 Confidence Intervalsfor
Percentiles

Many available computer programmes give estimates and confidence intervals (sometimes known as tolerance bounds) for percentiles,but their validity should be checked by comparison with reputable main-frame programmes. Approximate methods are available to calculate confidence intervals for percentiles, with the same theoretical basis as the intervals in 4.4, for those without access to avalid computer programme. Tables upon which the curves in Fig. 7 are based are available only for the first, fifth and tenth percentiles, but other tables are available in the literature. 451.1
Conjdence intervals for Weibull percentiles

The maximum likelihood estimate of thepth quantile, or ( 1OOp)th percentile ( j$ ), for the Gumbel distribution is: j$=C+hhln In _ 1 1-P H

iI

...(29)

Where u and B are the maximum likelihood parameter estimates. The confidence interval for y, is :
......(30)

yu(p,=iLGv,

......(31)

The maximum likelihood estimate of thepth quantile, or (lOOp)thpercentile ($r ), for the Weibull distribution 12

where yi @) andy, (p) are the lower and upper limits for the (100 p)th percentile. The Vi and I', factors for the relevant percentile are obtained from Fig. 7A and 7B.
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IS 11182 (Part 3/Set 2 ) : 1996 For the Gtmibel data in Table 2, for the fifth percentile, I', = 6.4 and VU=2.0, thus the 90 percent confidence limits amy (0.05) =4.1 kVandy,(O.O5) = 5.2 kV The maximum likelihood estimate of the fifth percentile is $,=4.95kV Thecon6denceboundsforthe pememiles, together with the confidence bounds for b, can be plotted on Gmnbel paper (Fig. 3B). `pable 1 BIW#&WIIData on Epoxy Specimens ~Clauses4.3.1.1,4.4.1and4.5.1.1)
Specimen No. i F(4) Breakdown Time r,

insulation on probability paper. However, this method can be misleading. This section outlines rigorous hypothesis test methods and presents an approximate simplified graphical procedme. 5.1 Rigorous Hypothesis Tests Various methods are available to test the hypothesis that one set of test results is `significantly' different from another. The following discussion is based on the Weibull distribution but, with the transformations in equations (4) and (5), it applies equally well to the It is difIicult to compare the distributions of two (or more) sets of data directly. Instead the following specific hypotheses, among others, are possible. Hp:P1 = I%!
H,:a,
Hp:xlp=

Percent 10 20 30 40
50

hours 15.3 30.3 48.5 89.4 90.4 105.7 144.9 _

1 2 3 4
5

......(32) ..... (33) ......(34)

6 7 8 9

60 70 80 90

=

CQ
XQ

where the subscript p refers to the pth percentile of the Weibull data, and the subscripts 1 and 2 to data sets No. 1 and No. 2. The tests of this hypotheses H, and HP often assume HP to be true. Rigorous means for proving or disproving the above hypothesis are usually based on `likelihood ratio' methods or &i-square (x2) approximations. Some tables are available for calculating the significance levels, but these tables do not cover all sample sizes or censoring. Computer programmes are available which are based on the likelihood ratio principle. Specific use of these programmes and the interpretation of the results will be discussed in the programme's documentation. For those without access to appropriate computer programmes, a simplitied `graphical' approach is presented. 5.2 Simplifwd Method to Compare Percentiles An approximate technique to determine if two data sets are different is to determine whether the confidence intervals at any percentile of the two distributions overlap, or not. If, for example, the confidence intervals at a specific percentile overlap, there is no indication that the two distributions are significantly different at that percentile. Conversely, if the confidence intervals do not overlap at the selected percentile, then it is likely that the two distributions are significantly different at that percentile. It is, however, ;m~nant to note that with small numbers of specimor;s the confidence intervals are, in general, relatively wide, so that the conclusions in such cases are rather weak. This type of test is independent of the shape parameter ofthe two distributions. Note, however, that ifthe shape parameters are different, a significant difference may exist at one percentile, but not at another. The conhdence intervals G.x the percentiles can be calculated as described in 4.5.1. 16

(Clauses
Specimen No.

Table 2 B-down Voltage of Oil 4.3.1.2,4.4.2,4.3.2.2 artd4.5.1.2) G(Yi)
Percent 9 18 27 36 45 54 64 73 82 91
Breakdown Yi(kV) 5.0 5.0 5.2 5.6 5.7 5.7 5.8 5.8 Voltage

2 3 4 5 6 7 8 9 10

5 COMPAIUSONS

A common situation involves testing two (or more ) insulation systems to determine which of the two is a `superior' insulation. Analysis of the results from such tests involves proving or disproving the hypothesis that there is a significant difference between the probability distributions of the data for the two systems. The `hypothesis test' is usually performed on the distribution parameters or at a given percentile, for example 63 percent, 50 percent or 3 1 percent p&abilities. For power systems insulation, it may be relevant to test for significant differences at low probabilities of failure, say 1 percent or 5 percent. The easiest approach to compare the results of tests on two types of insulation is to plot the test data for each
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IS 11182 ( Part 3/Set 2 ) : 1996 It is often useful to analyze comparison test data on Weibull probability paper. Plot the data from the two (or more ) tests on the same graph paper. As described in 4.5.1, plot the 90 percent tolerance bounds for each data set. Where the confidence intervals do not overlap, there is a high probability that the two data sets are significantly different. This procedure is illustrated with the test data in Table 3, which are also plotted in Fig. 8A. The 90 percent confidence interval for each data set is shown in Fig. SB. For percentiles greater than about 10 percent, the two intervals do not overlap and therefore there is a significant difference between the two insulations at high field stresses and at high probabilities of failure. Note that for low probabilities of failure the intervals overlap and thus a significant difference is not assured. In principle, more specimens need to be tested to determine if a significant difference exists at low probabilities of failure. Table 3 Breakdown Stresses (kV/mm) of Polyethylene from `Iivo Manufacturing Processes ( Clause 5.2 )
Failure No. Process 1 (Screened) 35 35 36 40 43 43 43 46 46 48 Process 2 (Unscreened) 39 4.5 49 49 53 53 53 53 55 55

and the laboratory test voltage ( or voltage stress ) for the specific insulation system tested, over the range of test conditions employed in the test programme. This laboratory performance relationship cannot be used direcly to establish intended performance, of the insulation system unrler service conditions in absolute terms. As is recognized in IS 11182 ( Part 1 ), it is desirable ultimately to be able to make evaluations, directly, in absolute terms. However, based on the present state of the art of electrical insulation technology, direct evaluation is not possible. Only comparative evaluations are possible at this time. If the estimated performance under service conditions of the identical insulation system that was subjected to laboratory testing can be evahmted in accordance with the principles outlined in IS 11182 ( Part 6 ) then a relationship between laboratory performance and estimated service performance can be established for this insulation system. If, in turn, the laboratory performance of a second insulation system is established by using an identical laboratory test programme, then the intended performance of the second insulation system under service conditions can be established by comparison with the relationship between the laboratory and estimated service performance of the first insulation system. The usefulness of these relationships is dependent on the validity and accuracy of the failure model chosen to represent them. 6.1 Choosing a Failure Model The selection of the most valild failure model for the test programme is equivalent to establishing the relationship between failure time and test voltage over the range of voltages covered by the several different constant voltage tests in the test programme. It is well known that, in general, elapsed test time in a voltage endurance test decreases as the test voltage stress increases. Although no individual model has been conclusively proven to be valid, two mathematical models have been most commonly selected to express this relationship quantitativeIy. Other models have also been used but will not be described here. 6.1.1 Inverse Power Model
The inverse power model is described by the equation: L = W-N .......(35)

1
2 3 4 5 6 7 8 9

10 11 12 13 14 15 16 17 18 19 20

48 48 48 48 48 51 51 51 51 57

51 57 57 57 61 64 64 65 67 68

6 INTERPREmTION ENDURANCEDA'IA

OF ELECTRIC&

The purpose of electrical endurance testing is to establish a relationship between the laboratory test life 19

where L is the time-to-breakdown at constant voltage I', and k-and N are constants. As the voltage increases, the time to breakdown decreases in inverse proportion to the voltage raised to the exponent N. A threshold value can also be incorporated.
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If the co-ordinates of time to breakdown and voltage (or voltage stress) fit a straight line when plotted on log-log graph paper, then the inverse power model and the equation (3 5) can be assumed to be valid. 6.1.2 Exponential
Model

as : L=Cexp[-(V-Vt)],V>V, ......(36)

The exponential model is described by an equation such

where C and Vtare constants. This model has a built-in threshold principle such that breakdown cannot occur at voltages ( or voltage stresses ) less than V,
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IS 11182 (Part 3/Set 2) : 1996 If the co-ordinates of time to breakdown and voltage ( or voltage stress ) do not fit a straight line when plotted on log-log paper, as described in 6.2, but do fit a straight line when the time to.breakdown is referred to a logarithmic scale and the voltage (or voltage stress ) is referred to a linear scale, then the exponential model and the equation (36) can be assumed to be valid. 6.2 Handling the Test Data The co-ordinates required for introduction into the selected failme model are time to breakdown and voltage ( or voltage stress ). The time to breakdown must be the characteristic or representative value derived from the calculations described in 4 for each level ofvoltage employed in the test programtne. It is beneficial to have the confidence intervals at the characteristic values for inclusion as well. 6.2.1 Pmgressive
Voltage Stress Tests

model has been assumed and calculated contidence intervals have been included. The voltage endurance data have been plotted on a log-linear graph in Fig. B3 and on a log-log graph in Fig. B4. 6.4 Interpretation The plot of time-to-breakdown against voltage (or voltage stress) and the exponent of the derived equation, as described in 6.3, become the figures of merit for the electrical endurance of the specific insulation system under the conditions of the test. 6.5 Limitations Any conclusions derived from statistical analysis of test data are rigorously valild only within the limits of the test conditions used to derive the data, and for the specific insulation systems included in the test pmgramme. Extrapolation beyond the limits ofthe conditions oftest is dangerous, and if a discontinuity exists in the relationship between failure rate and time, at different voltage stresses, extrapolation beyond that discontinuity is completely invalid and erroneous. In the absence of such discontinuity, the error of prediction ( even with a correct model ) becomes excessively large very quickly. 7 TESTP@OCEDURES AND REPORTS When insulation systems are tested, it is important for uniformity that the test procedures specify the following:

The test data resulting directly from the progressive voltage stress tests will be in the form of test voltage and rate of rise of voltage. For introduction into the selected failure model, these data must be transformed into coordinates of time to breakdown and voltage (or voltage stress). A number of procedures for making this mathematical transformation have appeared in the technical literature, and will not be described here. 6.2.2 Constant Voltage Stress Tests The test data resulting from constant stress will be in a suitable form for direct introduction in the selected faihue model. 6.3 Plotting the Data The co-ordinates derived in 6.2 are plotted in accordance with the selected failure model. Since the data provided by extreme value distributions at each test voltage are not symmetrical, it is not strictly valid to use a linear regression method to establish the equation for the straight line which best tits the relationship between those data on the voltage (or voltage stress ) versus time graph. An approximate procedure is to make a linear regression analysis of the logarithms of the time to breakdown data using any conventional arithmetic procedure or computer programme, and plot the antilogarithms of resulting line coordinates on the appropriate graph for the selected failure model. Confidence intervals at varying voltages (or voltage stresses) can-also be calculated and plotted. An example of the handling of a set of hypothetical data is given in Annex B. For this example, the inverse power

a) the Mrapolation model that should be used;

b) whether voltage or voltage stress should be
used;

c) whether voltage should be expressed in r.m.s or
peak values;

d) the voltage range or voltage stress range over
which the endurance values of insulation systems may be extrapolated, or the value to which the extrapolation should be limited;

e) the statistical treatment that should be used;
and

t) the specific values that should be used to
quantify performance. The report should include the plot of time to breakdown versus voltage (or voltage stress) and the numerical value of the exponent of the mathematical relationship, as well as the details specified by the equipment technical committee. These conditions and cosiderations may be different for various types of equipment and insulation systems.
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ANNEX A ( Clause 4.3.2.1 ) ( Informative )
BASIC PROGRAMME TO CALCULATE 8 et ^s 10 12 15 20 30 40 50 60 70 80 85 90 100 110 120 130 140 150 160 170 180 190 200 210 220 230 240 250 260 270 280 290 3(-)0 310 320 330 340 350 352 354 360 370 386 390 400 REM PROGRAMME CALCULATES MAXIMUM LIKELIHOOD REM OF WEIBULL PARAMETERS REM A FIRST GUESS FORBETA IS REQUIRED DIM TIME (29, LNTIME (29, A (3) PRINT `INPUT N-NO. OF SAMPLES TESTED' INPUT N PRINT `INPUT R-NO. OF SAMPLES FAILED' INPUT R PRINT `INPUT TS-TIME OR VOLTAGE TEST STOPPED' INPUT TS LTS = LOG (TS) PRINT `INPUT FIRST ESTIMATE FOR BETA INPUT BETA ITER=O c=o PRINT `INPUT THE FAILURE TIMES OR VOLTAGES' FORI=lTOR INPUT TIME (I) LNTIME (I) = LOG (TIME (I)) C=C+LNTIME(I) NEXT1 c=c/R FORJ= lT03 SUM=0 FORK= 1TOR SUM = SUM + ((TIME (K)"BETA)*(LNTIME(K)"(J-1))) NEXTK A(J) = SUM + (N-R)*(TS"BETA)*(LTS"(J-l)) NEXT-J QUOT = A (2) / A( 1) FPRIME = A(3) / A( 1) - (QUOT"2) + (( 1 ./BETA)"2) F=QUOT-(I./BETA)-C BETA =BETA - F/FPRIME ITER = ITER +l IF ITER>20 THEN 370 IF ABS (F)>O.OOO 1 THEN 200 ALPHA = (A( l)/R)"( 1./BETA) PRINT `ALPHA =`; ALPHA PRINT `BETA =`; BETA PRINT `ITERATIONS =`; TIER STOP PRJNT `DID NOT CONVERGE' PRINT `NEED A BETTER ESTIMATE OF BETA' STOP ESTIMATES

READY 22
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ANNEX

B

( Clause 6.3)
(Infomafive) MAXIMUM LIKELIHOOD

A WJCIDULLDISTRIBUTION TO MULTIPLE CENSORED DATA
FIT OF

Insulation system XYZ 8 kV/mm (r.m. s value) Test data Number of specimens tested = 10 Number of specimens failed = 10 Faihnetime 8 14 26 29 36 48 56 78 88 94 800 000 000 000 000 000 000 000 000 000

Suwival time Converged after 5 iterations output Alpha = 53582.97 Beta = 1.689 049 Insulation system XYZ 10 kV/mm (r.m.s. value) Test data Number of specimens tested = 10 Number of specimens failed = 10
*

Faihue time 150 280 490 800 1000 1250 1500 1600 2 100 3 000 Survival time Converged after 5 iterations output Alpha = 1340.608 Beta = 1.446667 23

IS 11182 ( Part 3/Set 2 ) : 1996 Insulation system XYZ 12 kV/mm (r.m.svalue) Test data Number of specimens tested = 10 Numberofspecimensfailed = 10 Failmetime 80 120 180 300 320 360 430 600 750 800 smvivaltime Converged after 5 iterations output Alpha=442.5018 EMa = 1.701543 Insulation system XYZ 14 kV/mm (r.m.svah@ Test data Number of specimen tested = 10 Numberof specimentiled = 10 Failuretime 4 6 7 12 17 19 24 29 33 40 Smvival time Converged after 5 iterations output Alpha = 2 1.422 23 Beta =1.680 513
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Insulation SystemXYZ8 kV/mm (n&s. Value) ( Co-ordinatesfir Figure B 1)
Sample NO.

status
System

VE hours

Rank

k

Hazard 100/k

i
I

n-(i-l)
10

Accumulated Hazard suM(1ooK)

1
2 3 4 5 6 7 8 9 10

XYZ
XYZ

F&d Failed Failed Failed Failed Failed Failed Failed Failed Failed

8 800 14 000 26 000 29 000 36 000 48 000 56 000 78 000 88 000 94 000

10.0 11.1 12.5 14.3 16.7 20.0 25.0 33.3 50.0 100.0

10.0 21.1 33.6 47.9 64.6 84.6 109.6 142.9 192.9 292.9

2 3 4 5 6 7 8 9 10

9 8 7 6 5 4 3 2 1

XYZ XYZ
XYZ XYZ XYZ XYZ XYZ XYZ

Insulation System XYZ lOkV/mm (r.m.s. Value) ( Co-ordinates for Figure B 1)
Sample Na Insulation System Status VJI hours 1 Rank i k Hazard Accumulated Hazard

n-(i-l) 10 9 8 7 6 5 4 3 2 1

100/k 10.0 11.1 12.5 14.3 16.7 20.0 25.0 33.3 50.0 100.0

SUM(iOO/K) 10.0 21.1 33.6 47.9 64.6 84.6 109.6 142.9 192.9 292.9

XYZ XYZ XYZ XYZ XYZ XYZ XYZ XYZ XYZ XYZ

Failed Failed Failed Failed Failed Failed Failed F&d Failed F&d

150 280 490 800 1 000 1 250 1 500 1600 2 100 3 000

1 2 3 4 5 6 7 8 9 10

2 3 4 5 6 7 8 9 10
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Insulation System XYZ 12 kV/mm ( r.m.s Value) ( Co-ordinatesfor
Sample No. Insulation System status VE

Figure B2 )
k Hazard Accumulated Hazard

Rank i

hours 1 2 3 4 5 6 7 8 9 10 XYZ XYZ XYZ XYZ XYZ XYZ XYZ XYZ XYZ XYZ Failed Failed Failed Failed Failed Failed F&d Failed Failed Failed 80 120 180 300 320 360 430 600 750 800

n-(i-l) 10 9 8 I 6 5 4 3 2 1

100/k 10.0 11.1 12.5 14.3 16.7 20.0 25.0 33.3 50.0 100.0

SUM (100 /K) 10.0 21.1 33.6 47.9 64.6 84.6 109.6 142.9 192.9 292.9

1 2 3 4 5 6 7 8 9 10

Insulation SystemXYZ 14 kV/mm (r.m.s. Value) ( Co-ordinates for Figure B2 )
Sample No. Insulation System Status VE Rank i k n-(i-l) Hazard Accumulated Hazard

hours 1 2 3 4 5 6 7 8 9 10 XYZ XYZ XYZ XYZ xiz XYZ XYZ XYZ XYZ XYZ Failed Failed Failed Failed Failed Failed Failed F&d Failed Failed 4 6 7 12 17 19 24 29 33 40

100/k 10.0 11.1 12.5 14.3 16.7 20.0 25.0 33.3 50.0 100.0

SUM (100 /`K) 10.0 21.1 33.6 47.9 64.6 84.6 109.6 142.9 192.9 292.9

1 2 3 4 5 6 7 8 9 10

10 9 8 7 6 5 4 3 2 1
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IS 11182 ( Part 3/Set 2 ) : 1996 Linear Regression of Logarithm of Time on Logarithm of Voltage Stress,y = II+ b x Insulation System : XYZ
( Co-ordinatesfir
Vdtage Stress UkV/mm Ig U(kv/mm)
X

Figure B2 )

Time

(t)

Igt
(s) Y 4.88 6.20 6.68 8.29 26.05

xy

xl

v'

hours

14 I2 10 8 c

1.15 1.08 1.oo 0.90 4.13

21 443 1 341 53 583

5.59 6.69 6.68 7.48 26.45

1.31 1.16 1.oo 0.82 4.29

23.80 38.47 44.67 68.65 175.59

Numberof observations:

n=4
r=

Degrees of freedom: cX=l()J
n CY

f=n-2=2

Average values:
Regression coefficients

' =20.17

y = T

= 6.51

a =y-bz

Dispersion of deviations ofy from regression line

s=

/E%pX.Y-~~~

=o.30

Standarderrorof u

s, = 2

= 0.15

Standard error of b

`b = de=

1.66

Forx=
Y= UkVlmm =

0.00 1.52 33.5

2.00 1.37 23.6

4.00 1.22 16.7

6.00 1.07 11.8

8.00 0.92 8.3

10.00 0.77 5.9

Corre~on

coefficient

r

==/&

=-o.g8
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